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Abstract 

Using  energy  conservation  and  causality  considerations,  the  completeness  of  scat¬ 
tering  states  in  established  for  plane  waves  impinging  on  an  irregular  interface.  Pro¬ 
vided  certain  limiting  operations  commute  with  differentiation,  it  is  shown  that  surface 
waves  need  not  be  explicitly  included  in  the  Weyl  representation  of  the  Green’s  func¬ 
tion  in  the  presence  of  a  rough  interface.  Rather  surface  waves  are  implicitly  included 
through  the  poles  of  the  scattering  amplitudes.  This  result  was  used  implicitly  in  a 
recently  developed  scheme  to  treat  scattering  in  a  duct  using  half-space  scattering 
amplitudes. 
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I  Introduction 


In  an  earlier  paper  the  method  of  smoothing  was  applied  to  wave  propagation  in  a  duct 
with  rough  boundaries.  The  idea  was  to  make  use  of  plane- wave  scattering  amplitudes  for 
each  of  the  surfaces  in  isolation,  i.e.,  scattering  amplitudes  for  each  surface  in  a  half-space. 

Fields  expressed  as  superpostions  of  scattering  states  associated  with  each  surface  were 
matched  across  a  fictitious  plane  passing  through  the  middle  of  the  duct.  Implicit  in  this 
scheme  is  the  notion  that  any  field  in  the  duct  below  (or  above)  the  fictious  plane  can  be 
represented  as  a  superposition  of  the  following  type: 

tf(r)  =  J  dQa(Q)  (exp(iQ  •  R  -  iqz)  +  J  dQ'[exp(*Q'  •  R  +  iYz)r(Q'|Q)]^  .  (1) 

Here  Q  represents  a  horizontal  wave  vector  and  q  =  \Ju2 /c1  —  Q2.  In  some  instances, 
however,  one  knows  there  are  surface  waves  which  decay  exponentially  from  the  interface. 

Are  these  included  in  this  formulation?  In  the  case  of  an  interface  between  two  fluids, 
how  does  one  incorporate  waves  propagating  up  from  below  the  interface?  The  aim  of  the 
present  work  is  to  examine  these  questions,  in  effect  to  show  under  what  circumstancee 
set  of  scattering  states  for  a  rough  fluid-fluid  interface  is  complete. 

As  indicated,  the  scattering  states  should  be  linear  in  the  scattering  amplitudes.  One 
learns  early  on  that  armed  with  a  complete  set  one  should  be  able  to  construct  a  Green’s 
function  for  the  wave  equation  bilinear  in  the  scattering  states: 

G(r.r')  «  J  AViJrJ^M/tk2  -  k").  (2) 

This  means  that  the  Green’s  function  will  contain  terms  bilinear  in  the  scattering  ampli¬ 
tudes.  On  the  other  hand,  the  Green’s  function  in  the  absence  of  a  boundary  is  expressible 
as  a  superposition  of  downgoing  plane-waves  below  the  source  level  (through  the  Weyl 
representation).  The  presence  of  a  boundary  can  be  accounted  for  by  converting  these 
downgoing  waves  into  a  sum  of  the  downgoing  wave  plus  reflected  upgoing  waves  which 
are  linear  in  the  reflection  coefficient  or  scattering  amplitude.  The  Green’s  function  in  the 
presence  of  a  boundary  is  thus  the  free-space  Green’s  function  plus  a  scattered  field  linear 
in  the  reflection  coefficients  [1].  How  does  one  reconcile  these  two  representations  of  the 
Green’s  function,  one  of  which  contains  terms  bilinear  in  the  scattering  amplitudes,  the 
other  of  which  contains  terms  at  most  linear  in  the  scattering  amplitudes?  A  hint  comes  0(les 
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from  the  optical  theorem  which  is  a  consequence  of  energy  conservation.  The  optical  the¬ 
orem  relates  the  imaginary  part  of  the  forward  scattering  amplitude  (linear)  to  the  total 
cross  section  (quadratic). 

In  this  paper  the  reconcilation  of  these  two  types  of  expressions  for  the  Green’s  function 
will  be  worked  out  explicitly  for  rough  interfaces  separating  two  fluids.  How  surface 
waves  and  evanescent  waves  come  in  will  be  shown  in  detail.  The  ideas  underlying  the 
reconciliation  are  energy  conservation,  causality,  and  time  reversal  invariance.  Of  course, 
these  notions  have  been  discussed  in  depth  for  quantum  mechanical  scattering  by  bounded 
objects  [2].  Here  the  scattering  surfaces  are  not  bounded  and  there  is  the  possibility  that 
the  form  the  wave  function  depends  upon  which  side  of  the  boundary  it  is  to  be  evaluated. 
These  features  make  the  interface  problem  slightly  different  than  the  quantum  problem. 
Additionally,  role  of  “closed”  channels  must  be  carefully  accounted  for  in  this  work.  In  the 
case  of  rough  interfaces,  (at  least  if  the  Rayleigh  hypothesis  holds)  the  wave  functions  can 
be  expressed  explicitly  in  terms  of  the  scattering  amplitudes  so  that  one  can  see  dearly 
how  the  connection  among  the  general  rules  plays  out. 

Earlier  work  treating  the  completeness  of  elastic  waves  in  a  solid  bounded  by  a  vacuum 
at  a  smooth  interface  was  done  by  Ezawa  [3].  He  showed  completeness  by  dose  examination 
of  the  special  form  of  the  reflection  coefficient  and  by  explidtly  considering  the  Rayleigh 
surface  wave  and  evanescent  waves.  Here,  the  motivation  is  to  show  how  general  prindples 
lead  to  completeness  relations. 

In  section  II  scattering  states  and  scattering  amplitudes  will  be  defined.  General  re¬ 
lations  such  as  reciprocity  will  be  reviewed,  particularly  the  generalized  optical  theorem. 
In  section  III  proof  of  completeness  will  be  presented  for  rough  fluid-fluid  interfaces.  In 
section  IV,  the  VVeyl  representation  of  the  Green’s  function  for  this  rough  interface  will  be 
discussed.  The  main  result  is  that  surface  waves  will  appear  explidtly  in  completess  rela¬ 
tions,  but  that  in  the  Weyl  representation  of  the  Green’s  function  with  an  interface,  poles 
in  the  scattering  amplitudes  account  for  surface  waves.  There  is  an  appendix  containing 
details  of  the  computations. 
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II  Scattering  states 

This  section  is  concerned  with  solutions  of  the  Helmholtz  equation 

(V1  +  (u»J/c!))V<r)  =  0  (3) 

in  two  half-spaces  separated  by  a  rough  interface.  On  the  upper  side  of  the  interface  the 
sound  speed  is  c+  and  the  ambient  density  is  p+.  Across  the  interface  there  is  a  sharp 
jump  discontinuity  in  these  parameters  so  that  below  the  interface  they  become  c~  and 
p~.  The  field  rp  represents  the  velocity  potential  for  particle  motions.  For  normal  fluids 
presure  and  normal  velocities  must  be  continuous  across  the  interface.  In  terms  of  the 
velocity  potential  these  conditions  are  expressed  by 

p+tp(r+)  =  P~rp(r~) 

n-'Vip(r+)  =  n-Vip(r~)  (4) 

(See  Ref.[4].)  In  addition  to  boundary  conditions  across  the  interface,  physical  solutions 
of  the  Helmholtz  equation  must  not  diverge  as  z  -*  ±oo,  however,  they  may  exhibit 
exponential  growth  in  limited  domains. 

Consider  scattering  states  associated  with  the  Helmholtz  equation  and  these  bound¬ 
ary  conditions.  For  each  frequency  u  there  are  solutions  (x+)  which  grow  out  of  plane 
waves  incident  from  above  the  interface  and  become  a  superpostion  of  plane  waves  leav¬ 
ing  the  interface  both  above  and  below  the  interface.  Likewise  there  are  such  solutions 
growing  out  of  plane  waves  incident  from  below  (x~).  Wave  vectors  here  have  horizontal 
components  denoted  by  uppercase  letters,  e.g.,  Q,K...  and  vertical  components  denoted 
by  corresponding  lower  case  letters  q,k....  (Similarly,  position  vectors  will  be  denoted 
by  upper  and  lower  case  letter,  e.g.  r  =  (R,  z)).  Because  the  plane  waves  must  sat¬ 
isfy  the  Helmholtz  equation,  vertical  and  horizontal  components  of  wave  vectors  are  not 
independent;  they  must  satisfy 

q±  =  y/w2/c±3  -  Q2  (5) 

above  or  below  the  interface.  The  coefficients  of  the  scattered  plane  waves  are  the  scatter¬ 
ing  amplitudes,  /2i,i(Q'|Q)  for  the  scattering  from  the  upper  fluid  (medium  1)  with  hori¬ 
zontal  wave  vector  Q  back  into  the  upper  fluid  with  horizontal  wave  vector  Q', 
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for  scattering  from  the  lower  fluid  (medium  2)  back  into  the  lower  fluid,  Ti^CQ'IQ)  for 
scattering  from  the  lower  fluid  into  the  upper  fluid  and  72,i(Q'|Q)  for  scattering  from  the 
upper  into  the  lower  fluid. 

Above  the  highest  point  of  the  interface  ( H+ )  and  below  the  lowest  point  of  the 
interface  ( H~ ),  scattering  states  can  be  written  as  follows: 


*o.a»  = 


exp(iQ  •  R  -  iq+z)  +  / dQ/exp(iQ/  •  R  +  *9+/*)i2u(Q'|Q)  if  z  >  H+ 
fdQ'exp(iQ'  •  R  -  if  z  <  H~ 


The  scattering  state  is  given  likewise  by 


fdQ'exp(iQ'  •  R  +  i?+,2)rlt2(Q/|Q)  if  z  > 

exp(iQ  •  R  +  iq~  z )  +  fdQ'  exp(t'Q'  •  R  -  *g”'z).R2)3(Q,|Q)  if  z  <  H~ 


In  these  equations  the  scattered  vertical  wavenumbers  q^1  are  given  by  Eq.5  with  Q 
replaced  by  Q'.  These  solutions  can  be  parameterized  by  u  rather  than  q±.  In  order  to 
produce  scattered  waves  which  travel  away  from  the  surface  or  which  decay  away  from 
the  surface,  if  q1  is  real,  it  must  have  the  same  sign  as  u,  and  if  it  is  pure  imaginary,  it 
must  be  positive  imaginary.  As  a  function  of  w,  q^  has  two  branch  points  at  iw/c*.  It 
will  be  convenient  to  take  branch  lines  extending  out  from  these  points  just  under  the 
real  axis.  Then  q  will  be  analytic  in  the  upper  half  complex  u  plane.  (See  Fig.  2).  This 
is  desirable  because  the  implicit  temporal  phase  is  exp(— iwt).  For  t  <  0  integrals  over 
u  can  be  closed  in  the  upper  half  plane.  Other  factors  should  be  analytic  there  so  that 
scattered  fields  vanish  before  the  incident  wave  arrives  at  the  scattering  surface.  (If  there 
are  surface  waves,  there  will.be  poles  in  the  upper  half  plane,  however). 

Recall  that  for  an  acoustic  field  with  pressure  p  and  particle  velocity  v  the  energy  flux 
is  [4] 

f  =  pv.  (8) 

For  fields  which  represent  time  harmonic  velocity  potentials, 

0  =  Re(^exp(-iwf))  (9) 


the  energy  flux  averaged  over  a  period  is 


fdt  =  {pu/  2)Im^*V^. 
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It  will  be  convenient  for  describing  some  of  the  formal  properties  of  the  scattering  ampli¬ 
tudes  to  use  plane  waves  normalized  to  have  unit  flux  in  the  vertical  direction: 

exp(tK  ■  R  ±  ikz)/yfpuk/ 2.  (11) 

It.l  Reciprocity 

Consider  horizontal  planes  above  and  below  the  rough  interface  separating  the  two 
fluids.  Let  4>a  and  ipb  be  any  two  solutions  of  the  Helmholtz  equation  satisfying  the 
b  -mndary  conditions,  Eq.4.  It  follows  from  Green’s  theorem,  the  Helmholtz  equation  and 
the  boundary  conditions  that 

[  dRp+(xl>adzipb-il>bdtil>a)=  [  dRp~ ( i>adz -  i>bdz rl>a ) .  (12) 

J Z^  J 

Actually  this  result  holds  even  if  there  are  a  number  of  fluid  or  elastic  layers  between  the 
planes  Z~  and  Z+ .  This  result  also  holds  when  the  continuity  of  normal  velocity  condition 
is  replaced  by 

n-VV>+  -  n-  Vi>~  =  -<7p+$+.  (13) 

For  flat  interfaces  this  last  boundary  condition  gives  rise  to  surface  waves  which  can  be 
seen  to  manifest  themselves  as  poles  of  the  reflection  coefficient. 

One  reciprocity  relation  follows  from  Eq.12  by  using  rf>a  =  Xq  and  V’i  =  X+j^: 

g+iM-Q|  -  K)  =  fc+.Ri,i(K|Q).  (14) 

By  using  V’i  =  XZ}£  one  ^n<^s 

p-*-r2tl(K|Q)  =  p+?+rlt2(-Q|  -  K).  (15) 

It  follows  from  using  tpa  =  Xq  and  =  x~j^  that 

g-i*2.2(- Ql  -  K)  =  *-A2i2(K|Q).  (16) 

n.2  Energy  conservation 

It  now  is  convenient  to  use  flux-normalized  plane  waves.  Define  incoming  plane  waves  by 

^±,Q,,±(r)  =  exP(*Q  •  R  T  *<l±z)/y/p*u 9*.  (17) 
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and  outgoing  plane  waves  by 


^±!g  ,,±(r)  =  exp(*Q  •  R  ±  iq±^)/yjp±^q±,  (18) 

where  the  signs  are  to  be  chosen  according  to  whether  z  is  above  or  below  the  interface. 

The  roots  in  the  normalization  are  to  be  understood  as  analytic  functions  with  branch 
lines  along  the  negative  real  axis.  Normally  the  sheet  that  will  be  used  will  have  non¬ 
negative  real  part.  Since  y/q*  =  y/q*  and  since  q  is  pure  imaginary  if  w/c  <  Q,  it  follows 
that  when  u  is  taken  to  be  an  independent  variable 

+  U(Q  -  «/c*)<Q^(r)  =  C-Q.,4(r)-.  (19) 

In  order  that  the  scattering  matrix,  to  be  defined  below,  reflect  reciprocity  as  simply  as 
in  Eq.26,  it  is  important  that  the  normalization  described  here  involve  q  rather  than  |g|. 
The  abbreviations 


0±(Q)  =  Ku/c^-Q)  (20) 

8*(Q)  =  0(Q-u/c±),  (21) 

will  be  used  in  the  next  section. 

Suppose  that  above  and  below  the  interface  for  a  limited  range  of  z,  arbitrary  super¬ 
positions  of  incoming  and  outgoing  flux-normalized  plane  waves  are  constructed.  Such  a 
field  above  the  interface  is 

0(r,u)  =  <.+(r)  =  j  dQ<.Q,,,(r)MQ)  +  ^ij^WoUQ).  (22) 

Below  the  interface  the  this  field  is  written 

V>(r,w)  =  V>-(r)  =  J  dQ<Qi,_(r)6.„(Q)  +  ,_(r)6014<(Q).  (23) 

Here  the  vertical  wave  numbers  are  determined  from  the  horizontal  wave  numbers  and  u; 
by  Eq.5.  It  is  not  claimed  yet  that  all  fields  can  be  represented  this  way;  that  would  be 
begging  the  question  of  completeness. 

If  these  states  obey  the  boundary  conditions  and  are  constructed  from  the  scattering 
states  x±,  it  follows  that  the  outgoing  coefficients  o^t  and  b^t  are  related  to  the  incoming 
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coefficients  a,n  and  6,„  through  the  scattering  amplitudes.  Because  normalized  waves  are 
now  being  used,  this  relationship  is 

(  ®.ut(K)  \  =  r  (  S,.i( K|Q)  5,,,(K|Q)  \  I  a.n(Q)  \ 

\  6„,(K)  J  •/  \  5j,,(K|Q)  S2l2(K|<J)  )  \  4,„(Q)  J  ’ 

where  the  scattering  matrix  5,  is  expressed  in  terms  of  the  scattering  amplitudes  by 

(  S!,.(K|Q)  Sj,j(K|Q)  \  /  ypji,.i(K|Q)  ^/pP^WQ)  \ 

V  S2.,(K|Q)  S2,2(K|Q)  )  ySFT!'l(KiQ)  J  ' 

In  terms  of  the  scattering  matrix  S  the  reciprocity  conditions  are  simply  stated  as 


S(K|Q)  =  S(-K|-Q)r. 


(26) 


Note  that  “anti-evanescent”  waves  may  be  included  in  these  superpositions  because  Q 
may  be  greater  than  u/c.  This  is  permitted  if  the  domain  in  which  the  superposition  is 
to  represent  a  physical  state  is  limited  in  z. 

To  see  how  energy  conservation  is  expressed  in  terms  of  5,  let  V’a  =  0*and  V’fc  =  V> 
in  Eq.12.  This  is  permissible  only  if  obeys  the  same  Helmholtz  equation  as  ,i.e., 
(cj/c*)2  must  be  real.  In  effect,  there  can  be  no  dissipation. 

By  taking  care  with  the  cases  in  which  vertical  wave  vectors  q±  are  purely  real  or 
purely  imaginary,  one  finds 

/  dQkn(Q)|2  +  /  dQlMQ)!2  -  2  /  dQIm(awt(Q)a,n(Qr)  = 

-/<9<|w|/c+  jQ<\w\/c~  JQ>|w|/c+ 

[  dQ|<w(Q)|2  +  [  <*Q| WQ)I2  -  2  /  dQM WQ)MQ)*)- 

JQ<  |w|/c+  JQ<|w|/c“  Jq>  M/c- 

(27) 

A  similar  result  for  elastic  waves  can  be  found  in  Ref.[[6]].  This  result  holds  for  arbitrary 
incident  wave  coefficients  a,„  and  6,n  and  hence  it  implies  constraints  on  the  scattering 
matrix.  To  separate  the  cases  of  propagating  and  non-propagating  waves  it  is  convenient 
to  define  the  matrix  operator  kernels  0  and  0  by 


0(Q|K)  =  | 

f  0(M/c+-Q)*(Q-K) 

0  V 

(28) 

^  0 

9(\u\/c--Q)6(Q-K)  J 

0(Q|K)  =  | 

(  &{Q  -  M/c+WQ  -  K) 

°  )• 

(29) 

^  0 

9(Q  -  \u\/c-)6(Q  -K)  ) 
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Now  because  Green’s  theorem  holds  for  arbitrary  incident  waves,  the  scattering  matrix 
must  satisfy 

0  =  S'OS  +  i{S'Q  -  0S).  (30) 

This  is  a  matrix  equation;  matrix  multiplication  is  implied  as  is  integration  over  interme¬ 
diate  wavevectors.  It  represents  a  generalized  optical  theorem  since  both  propagating  and 
non-propagating  waves  are  involved,  and  because  there  is  a  result  here  for  non-forward 
propagation.  A  result  of  this  type  was  used  by  Maystre  et  al  [8],  but  incident  “anti- 
evanescent”  waves  were  not  included  by  these  authors.  See  also  Ref.[7]  where  again  expo¬ 
nentially  growing  waves  were  not  considered. 

In  quantum  scattering  the  analogous  result  is  unitarity  of  the  S  matrix.  By  using  the 
reciprocity  relation,  Eq.26,  and  the  fact  that  the  matrix  operators  0  and  0  are  independent 
of  the  sign  of  Q,  the  horzontal  wavevector,  it  is  possible  to  show  the  following  alternative 
expression  of  energy  conservation: 

0  =  505f  +  i(Se  -  05f).  (31) 

In  a  space  in  which  all  waves  are  propagating  (all  channels  are  open)  these  relations  reduce 
to  the  unitarity  condition  in  quantum  scattering: 

1  =  55f  =  5f5  (32) 
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Ill  Completeness 


Using  just  the  scattering  amplitudes  it  will  not  be  possible  to  show  completeness  in  a  true 
sense.  This  is  because  the  representation  of  scattered  waves  as  a  superposition  of  outgoing 
plane  waves  weighted  by  the  scattering  amplitudes  might  not  converge  as  the  field  point 
is  moved  between  the  lowest  ( H~ )  and  highest  points  of  the  surface  (H+).  That  is  to 
say  that  the  Rayleigh  hypothesis  might  be  violated.  Nevertheless,  by  using  the  unitarity 
condition,  Eq.31,  it  is  possible  to  discuss  the  following  quantity  for  points  outside  the 
forbidden  region. 


A(r'r'>  =  (Sj5  Jdq  [r dq+xQ;*MxQ.^Tr  +  i°°‘ 7T&»xQ.,-MXQ,,-{rV 

(33) 

The  vertical  wavenumbers  associated  with  the  incident  wave  here  are  all  real.  Depending 
on  the  magnitude  of  the  outgoing  horizontal  wave  vectors,  the  scattered  vertical  wavevec- 
tors  may  be  either  real  or  pure  imaginary  since,  for  example, 


q+'  =  -  Q*  =  \/Q2  +  q(+)2  -  Q*. 


(34) 


In  the  remainder  of  this  section  it  will  be  shown  that,  in  fact, 


A(r,r')  =  *3(r-r'), 


(35) 


which  is  to  say  that  the  states  ase  complete  in  the  limited  sense  indicated  above, 
provided  there  are  no  surface  waves. 

The  quantity  A  is  evaluated  by  first  expressing  the  scattering  states  x  in  terms  of  flux- 
normalized  plane  waves  4>  and  the  scattering  matrix  S,  and  then  changing  from  integration 
over  vertical  wavenumbers  to  integration  over  the  frequency  u.  Since  the  incident  waves 
in  A  are  all  propagating,  horizontal  wavevectors  Q  are  restricted  to  Q  <  u/c^.  In  this 
fashion,  for  z  >  H+  and  z'  >  H+,  A  can  be  written  as 

f<x>  uip+ 

A(r,r)  =  /0  dW2(c+)2(2jr)3 


U 


dQ  [<Q.,+(r)  +  / dQV;“QV,(r)S,.,(Q'l<i:W)]  X 


Q<u//c+ 
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ii<„/c-iQV“QV'<r>5u<Q'lQ:">]  [/‘,Q>“Q"..*"(r')Sl-,(Q"IQi")]'}  (36) 

In  going  over  to  the  integration  over  frequency,  factors  of  c*-)2  were  canceled  by 

*-  -  m 

In  the  terms  quadratic  in  the  scattering  matrix  in  Eq.36,  one  can  recognize  the  (1,1)  part 
of  SOS*  that  appears  in  the  optical  theorem,  Eq.31,  if  the  integration  over  the  incident 
wavevectors  Q  is  performed  before  integration  over  Q#  and  Qw.  The  term  arising  from 
the  (1,1)  matrix  on  the  left  side  of  the  optical  theorem,  Eq.31  combines  with  the  terms  in 
Eq.36  to  give  5(r  -  r').  The  remaining  terms  are  linear  in  S,  coming  from  the  linear  terms 
in  Eq.36  and  the  linear  terms  on  t'.e  right  side  of  the  optical  theorem.  The  result  of  using 
the  optical  theorem  in  Eq.36  is  thus 

A(r.r')  =  Kr  -  r')  +  f  f  «f  J  «T 

{<r;;QV>)s»(Q'iQV)  [«+«")< ,QV>'>'  -  + 

[s+«5')<”QV.(r)  +  »*(0>3,V.(r)]  •  (38) 

Using  Eq.19  for  the  expressions  in  square  brackets  gives 

A(r,r')  =  <(  r  -  r')  +  j~  <^2(c“  ^  /rfQ'/rfQ" 

(39) 

Because  reciprocity  implies 

sJ^q'IQ'V)  =  5i>1(QwIQ,;«)  =  $u(-Q'l  -  QV).  (40) 

the  two  terms  in  braces  are  complex  conjugates. 

At  this  point  is  is  useful  to  switch  back  to  unnormalized  plane  waves  and  to  replace 
5i,i  by  R\t\.  Because  of  the  way  branch  cuts  have  been  taken, 

?(-u>)  =  -q*(w), 
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and  because  xq(r;u;)  and  X_q(r; -<*>)*  grow  out  of  the  same  incident  plane  wave,  one 
concludes  that 

*u( Q'IQ'»  =  RiA-Q'\  -  Q"; -«)*.  (41) 

Changing  signs  of  the  dummy  variables  Q'  and  Q"  for  the  second  term  in  braces  and 
letting  u>  — *  — u>  for  this  term  gives 

A(r,r')  =  i(r  -  r')  +  £  ^(^,5  / J *T 

exp(-iQ'  •  R  +  iQ"  •  R'  +  iq+'z  +  *?+,V)Rlli(Q,|Q,/;w)/?+,/.  (42) 

Anticipating  the  next  section,  one  might  identify  the  integral  above  to  be  the  time 
derivative  (because  of  the  factor  of  u)  of  the  scattered  field  arising  from  a  point  source 
at  r'  originating  at  time  t  —  0  and  evaluated  at  r  and  time  /  =  0.  This  is  what  one 
would  guess  from  the  Weyl  representation  of  the  free-space  Green’s  function.  Since  there 
should  normally  be  no  scattered  field  until  a  time  later  than  the  origin  of  the  incident 
field,  the  integral  in  Eq.42  should  vanish.  However,  if  there  are  surface  waves,  these  can 
exist  without  excitation  from  above.  In  that  case  the  “scattered”  field  might  not  vanish  at 
t  =  0.  The  existence  of  surface  waves  and  the  vanishing  of  the  integral  in  Eq.42  hinges  on 
the  analytic  properties  of  Rj,i  as  a  function  of  frequency,  u>.  Because  in  the  present  case 
z  and  z'  are  positive  and  the  vertical  wave  vectors  q  have  non-negative  imaginary  parts 
in  the  upper  half  complex  u>  plane,  one  should  be  able  to  add  a  large  semi-circle  in  the 
upper  half  complex  w  plane  to  the  contour  of  integration  in  Eq.42  and  evaluate  the  entire 
integral  using  the  residues  of  R,  assuming  that  it  is  possible  to  perform  the  integration 
over  oj  before  the  integration  over  Q'  and  Q". 

In  the  case  of  a  fiat  fluid-fluid  interface,  the  scattering  amplitude  is  given  by 

=  «!(Q  -  m 

The  reflection  coefficient  is  analytic  in  the  upper  half  v  plane  so  that  the  frequency  integral 
in  Eq.42  vanishes.  (Recall  that  branch  cuts  extend  outward  from  w+ie  =  ±c*Q  just  below 
the  real  axis.)  In  this  case,  the  integral  on  the  right  side  of  Eq.42  vanishes,  showing  that 
the  scattering  states  x*  are  complete. 

If,  on  the  other  hand,  the  boundary  condition  Eq.13  is  used  instead  of  continuity  of 
the  normal  derivatives,  and  if  c+  =  c~  and  p+  -  p~,  then  the  scattering  amplitude  for  a 
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flat  surface  is  given  by 


*.,i(QIQ»  =  <’(Q  -  (44) 

There  is  now  a  pole  in  the  upper  half  u-plane  at 

w  =  icjg2/*  -  Q2,  (45) 

provided  Q  <  g/2.  In  this  case,  the  integral  on  the  right  hand  side  of  Eq.42  becomes 

Wf  Ln exp(-'Q ■(R-R')-S(Z  +  <46> 

This  can  be  recognized  as  the  contribution  of  normalized  surface  waves 

XQ(r)  =  exp(-»Q-R-ff|z|/2)y|  (47) 

to  the  completeness  integrals. 

To  complete  the  demonstration  of  the  relationship  between  A(r,r/)  and  6(r  -  r'), 
A(r,r')  needs  to  be  evaluated  for  z  >  H+  and  z '  <  H~ .  This  is  done  in  the  Appendix. 
(Other  cases  of  z  and  z'  follow  from  these  two  iust  by  exchanging  the  indices  labeling  the 
two  media.)  The  result  is  that  for  z  >  H+  and  z  <  H~ ,  A  becomes 

£  exp(-iQ/  •  R  +  iQ"  •  R'  +  iq+'z  -  xq~"z')Tia{ Q'lQ"; u;)/q-w.  (48) 

This  integral  again  can  be  recognized  as  the  time  derivative  at  t  =  0  of  the  scattered  held 
evaluated  at  r,  above  the  interface  arising  from  a  point  source  located  below  the  interface 
at  r'.  If  the  scattering  amplitude  T  is  analytic  in  the  upper  half  u  plane,  the  integral 
vanishes.  Again  this  is  consistent  with  the  requirement  of  causality  when  there  are  not 
surface  waves. 

Thus  if  the  scattering  amplitudes  are  analytic  in  the  upper  half  plane, 

A(r,r')  =  £(r  -  r'), 

and  the  the  scattering  states  x*  are  complete.  If  so  the  scattering  states  here  are  the  nor¬ 
mal  modes  described  by  Tolstoy  [9].  Again  note  that  in  the  the  construction  of  A  incident 
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waves  are  all  propagating,  while  scattered  waves  may  be  propagating  or  evanescent.  Also, 
it  should  be  said  that  although  the  flat  surface  reflection  and  transmission  coefficients 
might  be  analytic  in  the  upper  half  u  plane,  this  is  no  guarantee  that  the  corresponding 
rough  interface  scattering  amplitudes  are  analytic  there.  Roughness  may  induce  surface 
waves,  at  least  for  the  mean  field  (See  [5]). 
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IV  Green’s  functions 


In  this  section,  the  results  of  the  last  section  will  be  used  to  construct  the  Green’s  function 
for  the  rough  interface  problem  at  least  outside  the  region  between  H~  and  H+.  Define 
A  a  similarly  to  A: 


1  r  f  da^c^)2 

Ac(r, r  ;  A)  =  J dQ  _  u(+)2)XQ,g+  ^r^Q,9+  (r  )  + 

/°°  dq~c <+)2  p+c^2  _ 

Jo  (X>-uI-)>)P-cMXQ«-VXQ.,-^ 


where 


w(±)2  =  c<±>2(Q2  +  9<±)2). 


(49) 

(50) 


The  Green’s  function,  G(r,r')  should  satisfy  the  same  boundary  conditions  across  the 
interface  as  the  scattering  states  and  in  addition, 


(V2  +  A 2/eW2)G(r,r';  A)  =  6( r  -  r'). 


(51) 


The  function  A g  will  be  examined  as  a  candidate  for  G.  It  is  dear  that  Ac  does  satisfy  the 
boundary  conditions  across  the  interface.  However,  it  is  also  clear  that  if  the  Helmholtz 
operator  can  be  passed  through  the  various  integrations  in  G,  then 

(V2  +  A2/c(±)2)AG(r,r/;  A)  =  A(r,r'; ) 

=<<r‘r')_£l  ssV(r’r';w)-  (52> 

The  function  0  is  an  abbreviation  for  the  double  integral  appearing  in  Eq.42.  When 
z  >  H+  and  z'  >  H+  then 

exp(-iQ'  •  R  +  tQ"  •  R'  +  iq+’z  +  *g+V)Ru(Q'|Qw;u,)/g+#'.  (53) 

When  z  >  H+  and  z'  <  H+,  ifr  is  read  out  from  Eq.48:  Rip  is  replaced  by  p+T\p/p~  and 
q+"  is  replaced  by  -q~+".  Note  that  reciprocity  implies  that 

V>(r|r';w)  =  V<r'|r;w)  (54) 
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Evaluation  of  Ac  can  be  carried  out  by  the  same  steps  used  to  evaluate  A,  keeping  track 
of  terms  quadratic  in  the  scattering  amplitudes  and  applying  the  optical  theorem,  Eq.  31. 
The  result  of  these  manipulations  is  that  A <3  becomes 


Aa(r,t;X)  =  (A,_ui) 

{jL*  (W'XP  liQ(R  '  R0  +  i,+(2 ' ^  /,+  -  i-)}  •  <») 

By  letting  A  — *  A  +  ie  and  adding  a  semi-circle  to  the  contour  in  the  upper  or  lower 
half  u  plane,  depending  on  the  sign  of  2  -  z\  the  first  integral  is  seen  to  be  the  Weyl 
representation  of  the  free-space  Green’s  function  with  frequency  A: 

G,„,(r|r';A)  =  /  (2T)4?,+  (A)  ***  liQ(R  ~  R>)  +  ~  •  <56> 

In  the  second  integral,  involving  the  contour  can  be  closed  in  the  upper  half  u  plane 
when  both  z  and  z'  are  positive.  There  are  two  kinds  of  contributions.  First  there  is  the 
residue  from  the  pole  at  u>  =  X  +  ie  and  second  there  are  the  contributions  from  the  branch 
cuts  or  poles  of  rp.  As  a  result,  A  a  can  be  written  as 

Ac(r|r';  A)  =  G/ree(r|r/;  A)  +  A)/2  -  jf  (A2  ^ 

The  contour  B  encloses  the  singularities  and  branch  cuts  of  if/  in  the  upper  half  plane 
but  not  A  +  ie.  When  acted  upon  by  the  Helmholtz  operator,  the  last  integral  yields  the 
integral  on  the  right  side  of  Eq.52.  Thus  if  this  integral  is  subtracted  from  A g»  the  result 
is  the  Green’s  function: 


G(r,  r';  A)  =  G/ree(r|r';  A)  +  V<r|r#;  A)/2  (58) 

One  can  check  easily  enough  to  make  sure  that  this  satisfies  proper  boundary  conditions 
across  the  interface,  if  this  construction  is  repeated  for  2/  below  the  surface,  and  if  the 
Rayleigh  hypothesis  holds.  The  point  of  this  exercise  is  that  the  poles  and  branch  cuts  of 
the  scattering  amplitudes  do  not  make  any  explicit  contribution  to  the  Green’s  function. 
By  adding  0  to  the  free  space  Green’s  function,  the  incoming  plane  waves  of  the  Weyl  rep¬ 
resentation  are  simply  replaced  by  the  scattering  states.  Surface  waves  are  automatically 
accounted  for  even  though  surface  waves  are  explicitly  required  for  completeness. 
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Appendix 


In  this  appendix,  A(r,  r')  is  evaluated  for  z  >  H+  and  z '  <  H~.  Some  the  the  details 
sketched  in  the  text  for  the  case  z!  >  H+  will  be  made  more  explicit  here.  For  example, 
in  terms  of  normalized  plane  waves  the  scattering  states  are  given  by 


xQ.,*(r)  = 


+  /■(QV£q-j,,Si,i(Q'IQ;w)]  z>h+ 

v/2fi/iQ'^v,SJll(Q'|Q;w)  z<H- 


XQ.,-(r)=  < 


(Al) 


V^/^QV0+“‘Q>i,+l(r)51,2(Q'|Q;u,)  z  >  H+ 

(A2) 

It  follows  from  the  definition  of  A,  Eq.  33  that 


A(ry)  =  (S?/ d^r iq*  Kq..*M  +  / X 

[/  <'Q'V~Q»^-.(r')ft..(Q"IQi«)]'  + 

(SF  [/ >< 

[<Q.,  -(«•)  +  /^'V“‘q.V„(r')52,2(Q"IQ;«)]'  (A3) 

Now  replace  the  variables  of  integration  by  u  using 

q±  =  yfu2l{c±)2  -  Q2  (A4) 

and  change  the  order  of  integration  so  that  the  integration  over  Q  <  cj/c*  is  first,  followed 
by  integration  over  Q'  and  Q"  and  then  u.  Multiplying  out  terms  in  S  gives 

yry>  =  (Spf 

i<^yQ/JQV“Q^«‘(t)s,^Q'IQ:"w-^-(rr+ 
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U  *“civ(r)Sl-J(Q'l<li“’)s»(<,l<,";"^“<JV"(r')'] }  •  (A5' 

The  integrals  over  Q  in  the  last  expression  are  equivalent  to  the  Q',  Q"  elements  of 

■5i,i0+(Q)51ti2  +  5ll20+(g)52>2. 

This  is  part  of  the  (1,2)  matrix  in  the  optical  theorem,  Eq.31  which  implies 

f  iShl(Q'\Q;u)Sl2(Q\Q";u)+f  Si^Q^S^QIQV)  = 

JQ<u)!c+  JQ<ut/c- 

-  i  [^(Q'lQV)©-^")  -  0+(g/)5,ti2(Q/|Qw;u;)] .  (A6) 

Substituting  this  result  into  the  last  expression  for  A  eliminates  the  terms  bilinear  in  S. 
The  remaining  terms  linear  in  5  can  be  written 

[^(flW'qV-W  +  «+(<3>“qV,(r)]  sUq'IQ";  ,.„w}  •  (A7) 

Usinq  Eq.19  to  replace  the  terms  in  square  brackets  gives 

f  J iV  /  *Q" 

’‘“'-q',,+.(r)'-fI.j(Q'IQ";")*~q"^-.(r')'} .  (as) 

The  reciprocity  relation  of  Eq.40  is  for  transmission  expressed  by 


Because 


it  follows  that 


■Su(QW»  =  ■s;.,(Q"IQ';“)  =  s;.,(-Q'l  -  Q";«). 


exp(— iQ  •  R  -  iq~(w)z)m  =  exp(*Q  •  R  +  iq~(-u)z), 


TiMlQ;  -u)  =  IT,2(-Q'|  -  Q;  -u;), 


(A10) 


(All) 
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and  therefore 

Si,2(Q'IQ";<*>)  =  Sf,a(-Q'|  -  Q";  -u).  (A12) 

With  this  result  the  two  terms  in  A  can  be  combined  to  give  a  single  integral  over  frequency 
from  -oo  to  +oo.  Expressing  the  result  in  unnormalized  plane  waves  and  finally  gives 

r*7$Rs?j**  i** 

exp(*(Q'  R  -  Q"  R')  +  iq-’z  -  tq-VjTUQ'lQVVq-'V).  (A13) 
This  is  Eq.48  of  the  text. 
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Figure  Captions 

Figure  1:  The  geometry  of  scattering  from  a  rough  interface.  The  rough  interface  sep¬ 
arating  fluids  of  density  p+  =  d+  and  sound  speed  c+  on  the  top  and  density  p~  =  d~ 
and  sound  speed  c~  on  the  bottom.  The  interface  is  bounded  between  parallel  planes  at 
z  =  H±.  Wave  can  be  incident  from  either  side  of  the  interface. 

Figure  2:  Contours  of  integration  in  the  complex  u;  plane.  The  dotted  half-lines  indicate 
the  branch  cuts  associated  with  q±.  The  dashed  line  and  the  real  axis  is  the  actual  contour 
of  integration  closed  by  a  semi-circle  at  infinity. 


20 


